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We isolated flux disorder effects on the transport at the critical point of the quantum magnetic field 
tuned Superconductor to Insulator transition (BSIT). The experiments employed films patterned 
into geometrically disordered hexagonal arrays. Spatial variations in the flux per unit cell, which 
grow in a perpendicular magnetic field, constitute flux disorder. The growth of flux disorder with 
magnetic field limited the number of BSITs exhibited by a single film due to flux matching effects. 

The critical metallic resistance at successive BSITs grew with flux disorder contrary to predictions 
of its universality. These results open the door for controlled studies of disorder effects on the 
universality class of an ubiquitous quantum phase transition. 


Transport phenomena near quantum critical points re¬ 
ceive ongoing scrutiny. Much attention comes from ef¬ 
forts to understand strange metal behavior in high 
mm and heavy fermion compounds mm- Others who 
are developing string theory based techniques to calcu¬ 
late many body properties of condensed matter in the 
strong coupling limit have focussed on quantum critical 
transport as well|5H2]. An interesting case is the super¬ 
fluid to insulator transition in charged two dimensional 
systems such as superconducting films. This Supercon¬ 
ductor to Insulator transition (SIT) appears in numerous 
thin film systems urn as a change from superconducting 
to insulating transport when a physical parameter such 
as thickness unmi] magnetic field [T^] or a gating field 
[IMTB] is tuned. At the critical value of the tuning pa¬ 
rameter, film resistances asymptote to a constant value 
in the zero frequency, low temperature limit. The limit¬ 
ing critical resistance assumes a value near the quantum 
of resistance for electron pairs ~ 4^2 = Rq [3 IS]- 
Despite the simplicity of this behavior, questions remain 
regarding its universality and the influence of disorder. 

Indeed, the critical resistance near a superconductor 
to insulator transition is at the heart of studies, both 
theoretical and experimental. Fisher m and collabora¬ 
tors m provided the original arguments supporting the 
universality of Re- They focussed on bosonic systems 
presuming that fermion degrees of freedom in the form 
of unpaired electrons play no role. They argued that 
i?c is likely to be constant in ordered systems within a 
universality class m especially in magnetic field. Dis¬ 
order appeared to lead to an increase in Re- Numerical 
simulations [siiin] also implied that quenched disorder 
modifies the SIT and Re- On the other hand, recent 
Quantum Monte Carlo simulations of ordered [7] and dis¬ 
ordered systems |2I] produced nearly equivalent values 
of Rc suggesting only a weak disorder dependence. Simi¬ 
larly, experiments produce conflicting results. For indium 
oxide films, Rc hovers around Rc — Rq [IH [22U24| for 
a large variation in microscopic disorder. Studies of ho¬ 
mogeneous amorphous bismuth films |25j . however, show 


that Rc grows with the normal state resistance, which 
is a similar measure of disorder. In both cases Rc as¬ 
sumes values that are about a factor of two higher than 
observed in the most ordered systems, micro-fabricated 
Josephson Junction Arrays (JJA) [53[ST]- Such an elab¬ 
orate picture reflects the variety of approaches and sys¬ 
tems employed to study Rc- It suggests that deriving a 
new method for controllably varying disorder and select¬ 
ing a thin film system with a bosonic SIT is crucial to 
making further experimental progress. 

Opportunities to carefully test quantum critical trans¬ 
port models of bosonic systems have arisen for thin film 
superconductors. There are clear indications that bo¬ 
son degrees of freedom dominate the SITs of a num¬ 
ber of them mm- Increasing resistivity or applying a 
magnetic field transforms superconducting transport into 
thermally activated insulating charge transport consis¬ 
tent with boson localization in for example InOa, |5S1 [53] , 
TiN [30], and nanostructured amorphous Bi films [3T1I32] . 
Magnetoresistance measurements showing Little-Parks 
oscillations [MIIMI indicate that the bosons. Cooper pairs 
of electrons, remain intact across the transition. Tunnel¬ 
ing experiments imply that the gap in the quasiparticle 
density of states persists into the insulating phase [5HI55] 
indicating that the fermionic degrees of freedom are not 
active at these SITs. Near the critical point, resistance 
data from the prototypical indium oxide film system show 
the expected scaling behavior [I3I131I3S]. Thus, recent 
work has made it possible to focus on bosonic phenomena 
using well chosen superconducting films. 

There are challenges to isolating the effects of disorder 
as it can take many forms that simultaneously influence 
film behavior at the SIT [53 EOl I57H55] . Models consider 
disorder as random variations in the electron potential 
[I7ll37l|38], or intersite coupling in lattice models [21] . 
or in physical parameters of grains in granular models 
[40] . Each of these can lead to qualitative accounts of 
SIT phenomena, such as the transition [8] , the emergence 
of granular structure in the Cooper pair distribution [41] , 
and the appearance of the giant peak in the magnetoresis- 


2 


tance of the insulating phase [55]. Distinguishing the in¬ 
fluences of each type of disorder, however, has been prob¬ 
lematic. The common disorder parameter, the normal 
state sheet resistance Rn, can reflect any of the above 
forms as it depends on carrier density, impurity poten¬ 
tial, and film morphology [8] . An additional confounding 
effect is that Coulomb interactions also grow with 
These repulsive interactions inhibit Cooper pair tunnel¬ 
ing between superconducting islands like those found in 
granular films [42] and microfabricated JJAs |43|. They 
strengthen as the interisland resistances increase through 
Rq to drive Cooper pair localization and thus, an SIT. 

We developed a method to study quantum critical 
transport of a bosonic system near the SIT that isolates 
the effects of one form of disorder from other forms and 
interaction effects. We achieved this control by creating 
systems with tunable flux disorder. The method employs 
thin films patterned with a disordered triangular array 
of holes. The multiply connected geometry of these so- 
called nano-honeycomb (NHC) films enables a single film 
to exhibit a series of magnetic field driven bosonic SITs 
[23| (see Fig. 1). The geometric disorder leads to vari¬ 
ations in the number of flux quanta per unit cell. This 
flux disorder grows with magnetic field so that suc¬ 
cessive SITs occur with increasing disorder. It limits the 
number of magnetic field tuned SITs that appear due to 
flux matching effects. More notably, rather than being 
universal, the Rc of these SITs increase with flux disorder 
from about 4 kD per square to plateau at about 6 kD/D. 
We discuss how this observation implies that geometric 
array disorder presumed to exist in microscopically disor¬ 
dered thin films HEum enhances Rc compared to ordered 
Josephson junction arrays. 

We used anodized aluminum oxide substrates with a 
nearly triangular array of holes |13| as a template for 
NHC films (Figs. la,e). All substrates had the same 
average hole spacing of 100 nm. The more strongly dis¬ 
ordered were produced by wrapping the aluminum with 
teflon tape to perturb the normally laminar flows that 
set up during anodization |4^. The geometric disorder 
is apparent in histograms of the unit cell areas for two 
typical NHC substrate (Figs. lb,f). The substrates were 
mounted to a dilution refrigerator and held at 8 K during 
film deposition. We studied arrays with fixed geometric 
disorder and varying normal state sheet resistance by de¬ 
positing a series of amorphous Bi films on a single sub¬ 
strate [33117]. Sheet resistances were measured at low 
frequencies using four probes. Transverse magnetic fields 
B were applied using a superconducting solenoid and are 
specified by the average number of flux quanta per unit 
cell in the array / = where A is the average unit 

cell area and $o is a flux quantum. 

Flux disorder results from variations in the geometry of 
the network of the templated NHC films. We characterize 
the network disorder by the fractional variation in the 
unit cell areas, 5a = A A/A where A A is the standard 
deviation calculated from gaussian fits to the unit cell 
area histograms (Figs. 1 b,f). In a magnetic field, there 


are variations in the local frustration, Sf = fSa that 
constitute the flux disorder [45]. This linear growth of 
Sf with magnetic field is presumed to dominate any field 
induced changes in other forms of disorder. In particular, 
previous work by our group m demonstrated that NHC 
films are comprised of an array of islands connected by 
weak links. It is supposed that randomness in weak link 
coupling and island size varies little with magnetic fields 
well below the estimated upper critical magnetic field [48] . 

The effects of flux disorder can be seen in a comparison 
of the low temperature magneto transport of two films 
with different geometric disorder but similar (Figs. 
1 c,g). In both cases, the magnetoresistance oscillates 
with a period of 1 between low values at integer / and 
high values at half-integer / ED). The oscillations decay 
more rapidly with / for the more disordered sample. In¬ 
vestigations of multiple substrates [H] indicate that the 
number of visible oscillations decreases from about 5 to 
1 as Sa increases from 0.03 to 0.14 and does not depend 
on R AT [49] . Thus, the data imply that oscillations persist 
only up to fields such that the flux disorder 6f « 0.3. The 
maximum number of oscillations observed in the most 
ordered arrays appears to be limited by the rise in the 
magneto-resistance that develops beyond 1 Tesla [3T] . 

Insight into the origin of the oscillations described 
above comes from prior experiments by Forrester et. al. 
EOl IS] on micro-fabricated Josephson junction arrays 
with geometrical disorder. Their magneto-resistances 
near the superconducting transition temperature Tc os¬ 
cillated with a period of 1 and decayed more rapidly in 
more disordered arrays. The oscillations completely dis¬ 
appeared above a critical field that was inversely pro¬ 
portional to the amount of geometrical disorder. The 
visibility of just 3 oscillations in Fig. Ig for 6a = 0.115 
is in rough accord with their results. In addition, they 
showed that the phenomenon results from Tc variations 
caused by oscillations of the average of the Josephson 
coupling energy Ej =< —J cos(0i — 9j — A^-) > between 
neighboring nodes in the array |50| . J is proportional to 
the amplitude of the superconducting order parameter 
on the nodes, 9i — 9j is the difference in the phases of the 
order parameter, and A^ = ^ A ■ dl is the line inte¬ 
gral of the magnetic vector potential between neighbor¬ 
ing nodes. Ej oscillates as A^- increases with magnetic 
field and the oscillation amplitude decays as the varia¬ 
tions in Aij grow. The resemblance of the oscillations 
presented in Fig. 1 with the prior results indicates that 
a similar modulation of Ej occurs in the NHC films. Its 
effect, however, appears more dramatic as near the SIT, 
Ej controls the strength of quantum phase fluctuations 
that potentially drive NHC films into an insulating state. 

Quantum superconductor to insulator transitions [52] 
are evident in the magneto-resistance oscillations in Figs. 
lc,g. They are identified by the crossing points in the 
Ra{B) traces taken at different temperatures. At each 
critical point (/c, Rc), the slope changes sign [44] as 
depicted by the i?n(T) in Figs. ld,h. Negative and pos¬ 
itive slopes correspond to films on the insulating and su- 
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FIG. 1. Comparison of NHC films on two disordered arrays, (a-d) show data from the more ordered array NHCl and (e-h) show 
data for the less ordered array NHC2. (a, e) Electron micrographs of the substrates. The unit cells, which were determined 
using a triangulation algorithm, are highlighted by yellow polygons. (b,f) Unit cell area distributions with Gaussian hts (red 
curves). The average unit cell area is j4 = 8 x 10^ nm^ for the two substrates, (c and g) Magnetoresistance oscillations 
(/ = BA/^o) at 130, 170, and 250 mK. (d and h) Film resistances as a function of temperature at zero held (/ = 0), near /c, 
and a.t f — 1/2. The hlms on NHCl and NHC2 have normal state resistances and thicknesses Rn = 17.9 kfl/D and d = 1.22 
nm and Rn = 19.2 kfl/D and d = 1.08 nm, respectively. 


perconducting sides of the SIT, respectively (cf. Ref. [8]). 
Changes in the slopes with field appear to be governed 
by smooth changes in an activation energy [Uj, which is 
consistent with the onset of boson localization. Also, 
the insulating film R[j{T)s show reentrant dips indica¬ 
tive of Cooper pairing fluctuations [33] prior to rising at 
low temperatures (Figs. ld,h). Seven crossing points are 
apparent in the more ordered sample while only three 
appear in the less ordered sample. 

It is also possible to observe a bosonic superconduc¬ 
tor to insulator transition at fields beyond the oscillation 
regime as shown in Fig. 2. The overlay of Ra{B) traces 
at different temperatures shows a crossing near 1.9 T. 
This critical field corresponds to /c « 8 or Sf « 0.9. 
Qualitatively, i?n(T)s at fixed magnetic fields develop 
tails at low temperatures that evolve into a flat depen¬ 
dence at the critical point with i?c ~ 6 kfl/D and finally 
an upturn with increasing field. These positive slopes 
dR\j/dT fit to an Arrhenius form, which implies an ac¬ 
tivation energy of these insulating films jUj. Also, these 
i?n(r)s show Cooper pairing reentrant dips. This evolu¬ 
tion resembles the SITs in the oscillation regime of thin¬ 
ner films (see Figs. Id, h) intimating that this high field 
SIT is also bosonic. The magnetic field induced reduction 
of Ej that drives this transition, however, must occur via 
a decrease in J caused by pairbreaking effects rather than 
an increase in frustration ^ A^j. 

The quantum critical resistances change systematically 
with flux disorder as illustrated in Fig. 3. The Rc ob¬ 
tained from films shown here and others [48] were deter- 



FIG. 2. A high field SIT. Film on NHC2 with Rn = 17 kH/D. 
(a) Isothermal magnetoresistance curves at 130, 170, and 250 
mK show a single crossing at Be = 1.9 T. Inset: Expanded low 
field region, (b) i?n(T) at discrete magnetic fields spanning 
Be- Vertical dashed lines indicate the isothermal slices in (a) 


mined from crossing points in continuous field sweeps as 
in Fig. 1 or by interpolating measurements of Ra{T) and 
at discrete fields to = 0 jUj. These methods 
yielded similar results when they could be compared. It 
is apparent in the inset that Rc increases with flux in the 
low flux limit for fixed geometric disorder. Linear fits to 
data from two individual films emphasize this monotonic 
rise. This behavior contrasts with ideal JJA’s for which 
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Rc is independent of flux. Plotting the Rc of a num¬ 
ber of films versus flux disorder 6f reveals that Rc rises 
with a slope of ^ 3 — 4.5 kfJ/D per unit of flux disorder 
to saturate near ^ for Sf > 0.3. Rc{5f 0) varies 
from 2.5 to 4 fcf]. Experiments on 1 set of films suggest 
that Rc(5f —>■ 0) depends on weak link coupling, which 
is unexpected for ideal JJA’s [45] . 

The flux disorder dependence of Rc in Fig. 3 provides 
an explanation for the difference in Rc measured in hlms 
and fabricated JJAs. Films showing bosonic SIT charac¬ 
teristics exhibit Rc « Rq with no clear dependence on 
Rn . Most notably, data on many different InOa, films in¬ 
dicate Rc « 5.8 kfl [l2l[22l [24| . Because these films lack 
clear geometrical structure, their BSITs presumably oc¬ 
cur in the high flux disorder limit. Experiments on JJAs, 
on the other hand, yielded 2.5 < Rc < 4.5 kfl |2^ and 
1.2 < < 2.45 kfl [27|. Clearly BSITs occur in the zero 

flux disorder limit in JJAs. Altogether these experiments 
point to Rc increasing with flux disorder (cf. Fig. 3). 



FIG. 3. Critical resistance as a function of flux disorder. 
Closed circles are high field SITs. Open symbols are SITs 
within the MR oscillation regime, with red circles from NHCl 
[Rn = 19.7 kfl/D) and blue symbols from NHC2: Diamonds 
{Rn = 20.3 kfl/D), stars {Rn = 19.2 kfl/D), open circles 
{Rn = 17.4 kD/D), squares {Rn = 16.7 kfl/D), solid circles 
{Rn = 16 kfl/D). Black symbols are from other substrates. 
Inset: Rc as a function of / for two of the films. Dashed lines 
are linear fits to the open circles of the same color. 


It is important to view the current results relative to 
previous experiments showing a relation between Rc and 
disorder as parameterized by the normal state sheet re¬ 
sistance Rn- Investigations of homogeneous MoGe hlms 
[53| and uniform amorphous Bi hlms |25j showed a strong 
positive correlation between Rc and Rn- This correlation 
could imply that Rc increases with disorder [25]. Alter¬ 
natively, it could be produced by fermion quasiparticles 
known to be present at these SITs |S3]. These fermions 
can provide a parallel dissipative channel that alters the 
transport in the vicinity of a nominally bosonic quantum 
critical point |SS|. Thus, the correlation between Rc and 


Rn could rehect changes in quasiparticle density rather 
than disorder [53]. We hasten to add that patterning 
these hlms with hole arrays can potentially yield insight 
into the interplay of bosonic and fermionic degrees of 
freedom at an SIT. It can reveal how dissipation effects 
on a bosonic SIT due to quasiparticles [55] evolve with 
increasing disorder. 

Although there have been many calculations of Rc at 
bosonic superconductor to insulator transitions includ¬ 
ing some that explicitly consider disorder effects [171112!“ 
milSSHSHI none appear to predict the behavior observed 
here. Uniquely, Kim and Stroud [5D] treated hux disor¬ 
der effects in Quantum Monte Carlo simulations of disor¬ 
dered square arrays of Josephson junctions. Their results 
predict that Rc decreases strongly with hux disorder from 
3.5Rq to 0.2iiQ for integer /, which is in sharp contrast 
to our results at non-integer / in Fig. 3. Likewise, a 
comparison [7] of a pair of Quantum Monte Carlo results 
[zun] employing the latest methods for extracting Rc 
from simulations suggests that Rc only weakly depends 
on disorder. In addition, the vast majority of predictions 
of Rc in the low disorder limit are higher than experimen¬ 
tal values by a factor of more than two PIZKIIlllIllssi- 
155] . This discrepancy could suggest that the universal¬ 
ity classes of the quantum rotor and Villain models em¬ 
ployed in some of the most advanced approaches mm 
do not match the experimental systems. Or, it could re¬ 
hect a need to further rehne methods for calculating the 
quantum critical conductivity in the zero frequency limit 
[31111113]. There is the further possibility that experi¬ 
ments have not probed the true quantum critical regime. 

In conclusion, we have described experiments reveal¬ 
ing that disorder inhuences quantnm critical transport 
of a bosonic system. We hnd that the critical resistance 
at the superconductor to insulator transition depends on 
hux disorder. This result invites more theoretical at¬ 
tention to its universality while illuminating a diherence 
between SITs in thin hlms and Josephson junction ar¬ 
rays. The results also highlight a lack of quantitative 
agreement between theoretical predictions and measure¬ 
ments of the critical resistance. The studies establish 
NHC hlms as uniquely suited for studying the effects of 
well dehned disorder on quantum critical transport and 
the universality class of a prominent quantum phase tran¬ 
sition. Further, they invite comparisons with other con¬ 
densed matter systems in which the interplay of quantnm 
criticality and disorder can arise such as doped strongly 
correlated electron systems [64] or cold atoms in random 
optical lattices [55] . 
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